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ABSTRACT 

The Galactic open cluster Westerlund 1 was found only a few years ago to be much more massive than previously thought, with 
evidence suggesting its mass to be in excess of ~10 5 M Q , in the range spanned by young extragalactic star clusters. Unlike those 
clusters its proximity makes spatially resolved studies of its stellar population feasible. It is therefore the ultimate template for a 
young, massive star cluster, permitting direct comparison of its properties with measurements of velocity dispersion and dynamical 
mass for spatially unresolved extragalactic clusters. 

To this end, we used the long slit near-infrared spectrograph VLT/ISAAC to observe the CO bandhead region near 2.29 pm, scanning 
the slit across the cluster centre during the integration. Spatially collapsing the spectra along the slit results in a single co-added 
spectrum of the cluster, comparable to what one would obtain in the extragalactic cluster context. This spectrum was analysed the same 
way as the spectra of almost point-like extragalactic clusters, using red superiant cluster members as velocity templates. We detected 
four red supergiants which are included in the integrated spectrum, and our measured velocity dispersion is 5.8km/s. Together with 
the cluster size of 0.86 pc, derived from archival near-infrared SOFI-NTT images, this yields a dynamical mass of 6.3xl0 4 M o . While 
this value is not to be considered the final word, there is at least so far no sign for rapid expansion or collapse. 
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1. Introduction 

For more than a dec ade now, following their discovery by 
IHoltzman et al.l dl992b . young massive star clusters have been 
the focus of intense studies, amongst other reasons because they 
represent excellent laboratories of star formation: their relative 
youth ensures that their present day mass function still closely 
resembles the initial mass function (IMF), their high mass in 
principle allows studies of IMF variations (stochastic sampling 
of the IMF plays a minor role in clusters with masses in excess 
of approximately 1O 5 M ), and last, but not least, they are abun- 
dant in different types of galaxies where they probe a variety of 
environments. They suffer only one major drawback, the large 
distance, which means that they are slightly resolved at most, 
and typically also fairly faint. 

The ideal candidates for star formation studies would show 
the advantages of the extragalactic clusters, but sit in the Galaxy 
and consequently allow spatially resolved analysis of individual 
stars in the clusters. Westerlund 1 (Wd 1) is potentially such an 
ideal target, the first one discovered in the Milky Way (at a dis- 
tance of 5 kpc, where 1" corresponds to 0.025 pc). Even though 
it was already detected in the 1960's ( Westerlund! 1 1 96 lb . it took 
a few decades to find out that the cluster light is highly extin- 
guished by large quantities of dust, and that it possibly has a 
fairly high mass of about 1O 5 M . 

We observed it as a template for the almost unresolved ex- 
tragalactic star clusters. For most of those which were analysed 
spectroscopically, it is plausible that a large fraction of the clus- 
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ter mass is contributed by low-mass stars (e.g.. IMengel et al.l 
(l2002l) ; lLarsen & Richtler!(l2004 : lLar"sen et alj 32004b ). and that 
the clusters are compact enough to survive for at least sev- 
eral Gyrs. On the other hand, age determinations of individual 
clusters in a cluster population suggest that around 90% of all 
young massive extragalactic st ar clusters, irrespective of mass , 
are disrupted within 10 Mvrs jFal 12004 IMengel et all 120051: 
Bastian et al., 2005; Whitmore, 2006). Of course, it is possible 
that the clusters targetted for spectroscopy have already under- 
gone a selection effect, because they are all older than 6-7 Myr 
and therefore clearly not amongst the first ones to disrupt. Other 
types of investigation are necessary to test this hypothesis. 

But the origin of this discrepancy could also lie in false as- 
sumptions concerning the dynamical state of the young extra- 
galactic star clusters - for example, disruption or collapse could 
lead to an increase in velocity dispersion which is unrelated to 
cluster mass. Therefore we wanted to test the approach used for 
extragalactic clusters on the only spatially resolved young star 
cluster which has a mass around 1O 5 M . 

2. Observations and data reduction 

Observations were conducted on VLT-ANTU on the nights of 
March 1 1 and 12, 2006 and consisted of scanning the slit across 
a region while integrating. The ISAAC slit has a length of 2', 
and we selected a slit width of 073 to obtain the highest possible 
specttal resolution. This is necessary for determining a stellar ve- 
locity dispersion below lOkm/s, which we expected for the clus- 
ter. The scans were conducted by acquiring a star in the centre 
of th e cluster (W33 - thr oughout the text, we use the numbering 
from Wes terlundl (Q~987)), offsetting 2' East, and then scanning 
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Fig. 1. Top: 2MASS K-band image of the centre of Westerlund 1 . 
Indicated is the region over which the slit was scanned. Three 
scans were performed, each starting on the eastern edge and 
scanning the slit across the field of view during the five minute 
integration. The brightest stars are saturated in the image. Small 
black full circles are late-type stars, small black broken circles 
are the NIR-brightest early type stars, and the large broken red 
(grey) circles indicate the half-light radius (including uncertain- 
ties). Bottom: Mosaic of the three spatial scans of Westerlund 
1, integrated over 4' perpendicular to the slit. Increasing wave- 
length from bottom to top. Sky subtracted, flatfielded and dis- 
tortion corrected, aligned with image above. This 2d frame was 
collapsed along the x-axis in order to obtain the Id spectrum. 



westward over 4', for a total integration time of 300 s. This was 
repeated three times, and two more scans were conducted by off- 
setting 108" North and South, respectively, from the acquisition 
position. Figure Q] shows the region covered and the resulting 
2-dimensional spectrum. 

The central wavelength was set to be 2.31 /mi, and the actual 
wavelength range covered is 2.249 - 2.373 /mi at a resolution of 
R*9000. 

Reduction of the data was performed using the IRAF pack- 
a gfl For the slit scan, it included sky subtraction (a scan across 
an empty sky field 10' away was performed with identical inte- 
gration time of 300 s directly after the cluster scan), bad pixel 
correction (but there was a blotch of bad pixels (dust) in the top 
left quadrant which appeared during the night which neither flat- 
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fielded out, nor was properly removed by median interpolation), 
and flatfielding. 

Wavelength calibration combined the inf ormation of the 
XeAr lamp spectra and the night sky lines (Rous selot et all 
2000). Spatial distortion was easily corrected due to the many 
star spectra along the slit. The three scans (centre, offset North 
and offset South) were mosaicked into a combined scan before 
integrating spatially into a single cluster spectum. 

Individual spectra of the four near-infrared brightest stars 
(indicated with small solid circles in Figure [TJ, which show 
prominent CO absorption, were extracted as template stars for 
the determination of the cluster velocity dispersion. 

Since solar type stars were observed as telluric stan- 
dards (Hip077718 (G2V, 111^=5.894) and Hip087370 (G3V, 
rns:=5.609)), we divided them by an adequately smoothed high- 
resolution solar spectrum before applying the correction of tel- 
luric features to the Westerlund 1 spectrum. The features near 
the Mgl line around 2.2814/zm appear to be a residual from this 
step which we could not eliminate. 

Finally, we normalized all spectra which were used for the 
analysis by dividing by the continuum level determined from a 
small window (roughly 50 pixels) just shortward of the CO band- 
head. A rough flux calibration was obtained via a two-step ap- 
pro ach using the archival SOFI/NTT narrow band data obtained 
for lCrowtheretal.1 d2006l) . programme ID 073.D-0321, and the 
relative intensities from our scan. Even though we used a com- 
bination of NB 2.09//m and NB 2.17/mi, we could reproduce 
the magnitudes published by Crowth er et al.l J2006) to within 
0. lmag, and determined approximate K-band magnitudes for the 
moderately bright cluster members not listed in their publication. 
For the brightest cluster members, we could not use this image, 
because they were saturated. In those cases we used relative in- 
tensities to fainter stars common to the scan and the image (and 
isolated enough such that spectra of different clusters were not 
projected on top of each other). The relative fluxes of the stars 
remained quite constant throughout the short scanning time of 5 
minutes, because the seeing varied only slightly. The three indi- 
vidual scans of the central region showed variations in the rela- 
tive stellar intensities of up to 20%, more typically around 10%. 

3. Analysis and Results 

The resulting spectra are shown in Fig. [3] We additionally in- 
dicate the region which was used to determine the equivalent 
wi dth of the 2.29/«n CO bandhead according to the definition 
of lOriglia & Olival (2000). From this equivalent width Wco2.29, 
in combination with approximate absolute magnitudes, we deter- 
mined the rough spectral type of the individual stars as explained 
in Fig. [2] 

For the K-band magnitudes , we u sed a distance modulus 
of 13.5, my from IWesterlundl dl987l) . and apply an extinc- 
tion correction of 0.7 mag in K-band, which is a lower limit 
(ICrowtheret all 120061) . 

Fig. |2] shows our measurements of Wco2.29 in comparison 
with published values as a funct ion of T e f f . The lower par t of the 
figure is based on a figure from lForster-Schreiberl (120001) which 
includes measured values for Wco2.29 for giants and RSGs from 
her own measurements and two other publications (see refer- 
ences in that work). The loci populated by giants and SGs are 
indicated as shaded areas, and our measured values for both, 
Wco2.29 and are represented as horizontal bars which also 
indicate the measurement uncertainties. The range of effective 
temperature for which a given measurement overlaps with mea- 
sured values for stars is indicated by the shorter (labelled) hor- 
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izontal bars in the top of the lower part of the figure. In the top 
part, we applied essentially the same technique. F or theoretical 
Mk, we used V-K colours from | Koornneefl dl983l) . and theoret- 
ical My from lSchmidt-Kalerl d!982l) ). For three stars, both diag- 
nostics lead to an agreement in spectral type / T e fj. For the two 
brigh test stars, the mismatc h is caused by the lack of data points 
in lForster-Schreiberl d2000t) below 3000 K. Therefore, the table 
lists the spectral type resulting from Nig- . 

The four brightest stars are assigned spectral types ranging 
between M3I and M6I (foreground objects are ruled out also 
because of the small spread in radial velocity). In order to de- 
termine if W263 is a member of Wd 1, and hence valid for in- 
clusion in our velocity dispersion estimate, we determined its 
spectral type and radial velocity. It is much fainter than the red 
supergiants, and we considered Wco2.29 and V- and K-band pho- 
tometry to decide if it is a cluster member or not. The equivalent 
width would be consistent with a late K-type supergiant or a late 
M giant. But the absolute luminosity is much too faint, an M5III 
star approximately at the distance of Wd 1 is likely, or a similar 
spectral type if it is somewhat closer or more distant. Moreover, 
its radial velocity is very different from that of the other stars; 
the velocity difference of around 60km/s would be several times 
the escape velocity of the cluster. Hence, we have rejected W263 
from our dynamical analysis. 

We compared Wco2.29 = 8.7 A of the integrated spectrum 
with predictions from Starburst99 (Vazquez & Leitherer, 2005), 
using an instantaneous burst at solar metallicity with a Kroupa 
IMF (a = 1.3 between .1 and 0.5 M Q , a - 2.3 between 0.5 and 
100 M , iKroupal (|2002|) ). The modelled W C0 2.29 reaches 8.7A 
at an age of 6.3 Myr. It is difficult to quote uncertainties for this 
age, because the main uncertainty is not caused by measurement 
uncertainties in the equivalent width, but rather in the models 
(e.g. treatment of the post-main-sequence phase), the assump- 
tions (instantaneous burst), or simply statistics. Therefore, while 
the uncertainty in the equivalent width would lead to an age un- 
certainty of only around 0.2 Myr, the real age may easily be off 
by one or two Myr. Some of the spectral properties of the clus- 
ter as a whole and of individually extracted stars are listed in 
Table 1. 

The integrated spectrum was used together with the indi- 
vidual spectra to determine the velocity dispersion in the same 
manner as with spectra of extragalactic star clusters: broadening 
the continuum-diluted red- or blueshifted template spectra^ by 
convolution with a Gaussian with a given velocity width until 
a best fit (determined through a x 2 minimization) of the cluster 
spectrum is ob tained. From th e estimated total mass of the clus- 
ter (10 5 Mn. [Clark et all d2005l) ). we expected a velocity disper- 
sion of around 9 km/s, still quite comfortably above the limit of 
around 7 km/s to which we can determine a velocity dispersion 
in that spectral ran ge if the signal-to-n oise ratio is high enough 
(see Fig. 4, top, in lMengel efatl d2002l) ). 

However, our measured velocity dispersion is much lower 
than expected. For different template stars (which were of course 
the individual stars extracted from the scan) and different wave- 
length ranges, our fitting routine returned results between 1 and 
6 km/s. As there was no trend in our determined value with ei- 
ther spectral type of the template or spectral region considered, 
it was obvious that this was the random scatter expected for a 
velocity dispersion below the resolution limit, and we can only 
deduce that the velocity dispersion is below 7 km/s. 



2 Velocity dispersions from broadening of absorption features can 
only be measured for clusters dominiated by a stellar type with pho- 
tospheric absorption and without fast rotation 
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Table 1. Properties of the integrated cluster spectrum and of the 
five late-type stars included in that s pectrum. The columns list: 
Designation as in IWesterlundl (fl987h . W C0 2 29 according to the 
definition in Origl ia & Olival d2000h . the apparent and absolute 
K-band magnitudes (where the uncertainty is around 0.3 mag 
for the first, and 0.4 mag for the latter) and the spectral type of 
the stars determined from Wco2.29 and the absolute magnitude. 
For the cluster, the last column lists the cluster age determined 
from Wco?.?.9 in compariso n to predictions from Starburst99 
d Vazquez & Leithererll2005l) . The apparent magnitude listed for 
the integrated spectrum is only given in order to be able to com- 
pute the relative contributions of the individual stars to the to- 
tal spectrum from it. W26, the NIR brightest object, contributes 
20%, and all the late type stars together around 40%. The four 
brightest stars have magnitudes identifying them as red super- 
giants at the distance of Wd 1 . 
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Fig. 2. Determination of the spectral types of the observed stars. 
Our measured values for Wco2.29 and Ma: are compared to pub- 
lished measurements as a function of T„ff (|Forster-Schreibe"rl 
2000] iKoornneefl 119831: ISchmidt-KaleR Il982h~ For details, see 
text. 
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Fig. 4. Spectra of the brightest red supergiants and the red giant 
not corrected for radial velocity. The lines (bold: CO 2-0 band- 
head and CO(2-0) R13 rest wavelengths, dashed: centerd on two 
arbitrary absorption lines in the lowest red supergiant spectrum) 
are supposed to guide the eye, showing that a relative velocity 
between the stars can easily be determined. Below the names, 
heliocentric ally corrected radial velocities are given. 

What is possible to determine, though with relatively high 
statistically uncertainty, is the rms radial velocity dispersion of 
the four red supergiants. The rms velocity dispersion determined 
from the - unfortunately only four - individual stellar veloci- 
ties is 5.8+2.1 km/s. Weighting by the luminosity lowers this 
value to 5.0+1.7km/s. Errors on both of these values were de- 
termined through Monte Carlo simulations, using the uncertain- 
ties in the individual velocity determinations. Formally, the lat- 
ter value would be the one to compare with values obtained from 
observations of extragalactic clusters. But since this comparison 
is inappropriate in this case, due to the fact that Westerlund 1, be- 
cause of its young age, contains far fewer red supergiant (RSG) 
stars than all extragalactic clusters studied so far, we use the dy- 
namically more meaningful unweighted rms. Apparent magni- 
tudes of the stars are given in Tab. 1, and the heliocentrically 
corrected radial velocities in Fig. |4] 

The other parameter we require for the determination of the 
dynamical mass is the projected half-light radius. We used the 
NIR image we created from the archival NTT/SOFI data to ob- 
tain a radius which contains half of the total light (in projection) 
to be r/, p =0'.'86+0'.'14. The uncertainty comes mostly through as- 
sumption of different locations of the weakly constrained cluster 
centre. 

The resulting dynamical mass, for a cluster in Virial equi- 
librium, is (assumin g 77 = 9.3, for details on the formula see 
iMengelet all (120021) 1 

M dyn = ■qa 2 r hp IG = 2161 x 5.8 2 x O.86M = 6.3^ x 1O 4 M 

Within the e rrors this is consi stent with the mass expected from 

photometry dClark et all 120051) . 

From stellar number counts, we also (like lClarket ^ (12005)) 
arrive at an expected total cluster mass of approximately 1O 5 M . 
We see 4 red supergiants, a nd take the number of Wolf-Rayet 
stars (WN: 16, WC: 8) from lCrowther et al l j2006h , and of blue 
supergiants (BSG; > 25) from IClark et al.l d2005l) . The relative 
numbers of these stellar populati ons can be used tog e ther w ith 
evolutionary synthesis models (Vazquez & Leitheren, 120051) to 
age date the cluster. The model we assumed was the same we 
described earlier in this section. From this, according to Fig. 
[5] we determine the best fitting age to be 5 Myrs. The abso- 
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Fig. 5. For a cluster model generated with Starburst99, v.5.0 
(Vazquez & Leitherer, 2005), we assumed solar metallicity, an 
instantaneous burst of a 10 6 M Q cluster, and a Kroupa IMF. For 
cluster ages between 4 and 5 Myrs, we determined the expected 
star ratios from the output HR diagram in time steps of 0.1 
Myrs (after 5 Myrs in steps of 1 Myr). These theoretical values 
were, for each time step, divided by the observed ratios deter- 
mined from these absolute star counts: RSG=4, BSG>25, W- 
R-stars: 24, This means that for each individual diagnostic line, 
the model age is where this normalized ratio is 1 . For all three 
diagnostic lines this occurs at 5 Myrs. 

lute stellar numbers then indeed correspond to a 10 5 M Q cluster. 
Uncertainties on cluster mass from this method are difficult to 
determine, because of the sudden appearance (or disappearance) 
of post-main-sequence stars, RSG and Wolf-Rayet stars being 
almost mutually exclusive. Nevertheless, if the cluster age was 
wrong by one Myr, the total cluster mass could be as small as 
5xlO 4 M . One puzzling feature remai ns that sup e rgia nt stars of 
spectral type B3Ia were reported in IClark et afl d2005l) . which 
according to the model only appear at a cluster age of around 
12 Myrs. This could indicate that the assumption of an instanta- 
neous burst is invalid. 

4. Conclusions and future work 

1 . We identify four late-type supergiant stars (M3TM6I) which 
in total contribute 40% of the light in K-band in the scanned 
region. The equivalent width of the CO bandhead of the inte- 
grated cluster spectrum is 8. 7 A, which translates to a model 
age of 6.3 Myrs, somewhat older than the 5 Myrs determined 
from the stellar content. This discrepancy could arise from 
uncertainties in the models (the same model is used for both 
methods of age dating), or, certainly in part, because some 
of the stars which would contribute to the diluting contin- 
uum are not included in the spectrum because they are either 
outside the scanned region, or extinguished by dust. 

2. The velocity dispersion as determined from the integrated 
spectrum is below 7 km/s; the rms velocity dispersion from 
individual stars was found to be 5.8+2. 1 km/s. The resulting 
dynamical mass is 6. 3+37 x 10 4 M . Of course, these results 
cannot be considered final, because there are at least three as- 
pects which could influence the velocity dispersion. Firstly, 
of course the uncertainties are rather large for an rms value 
determined from only four stars, and we need more stars for 
a more reliable estimate. Secondly, mass segregation is not 
ruled out in this cluster, since the main sequence has not yet 
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Fig. 3. Specta of integrated cluster and the near-infrared brigthest individual stars. Atmospheric absorption (corrected for, but over- 
compensation a nd sub-optimal sky subtraction causes a few artefacts at long wavelengths) is shown at the bottom. Line identifica- 
tions taken from lHinkle et al.l (Il995l) . 



been detected. In that case, even measuring more stars would 
yield a velocity dispersion which is too low for the total clus- 
ter mass, because the mass segregated high-mass stars have a 
smaller-than-average velocity, and also velocity dispersion. 
Furthermore, these most luminous cluster members would 
cause an artificially low estimate of the cluster size. In that 
case we would have a clear demonstration of the perils of di- 
rectly interpreting stellar velocity dispersion measurements 
as reliable mass indicators in more distant clusters. And 
thirdly, this cluster is at the lower mass limit where observ- 
ables become relatively unaff ected by statistical fluctuations 
in the sampling of the IMF (ICervino & Valls-Gabaudl l2003t 
ICervino & Luridiana} 120061) . which is why we observed it. 
Nevertheless, its probable mass, and hence the number of 
high mass stars, is only 1/10 that of some of the high mass 
extragalactic star clusters. 

3. There is not yet a sign for rapid expansion or collapse, which 
would have shown up as an unexpectedly large velocity dis- 
persion. The assumption of Virial equilibrium thus appears 
to not be strongly violated. 

4. As this observational technique and analysis works well, fur- 
ther work at a wavelength more sensitive to the far more nu- 
merous hotter stars in this cluster would greatly increase the 
significance of our preliminary findings. 
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